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ABSTRACT 



Gravitational microlensing of quasars by stars in external galaxies can introduce fluctuations 
in the centroid of the "point-like" macro-images. The induced shifts are extremely small, on 
micro-arcsecond scales, below the limits of current optical observations. However, such shifts 
will become measurable with the proposed "Space Interferometry" mission, due to fly in 2005. 

The degree of the centroid shifts and their application as probes of both quasar structure and 
the stellar mass function in the lensing galaxy are discussed. 



Subject headings: Gravitational Microlensing; Quasar Structure; Interferometry 
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1. 



Introduction 



Over the last decade gravitational microlensing has shown itself to be not only a probe of the mass 
structure in our galactic halo (Alcock et al. 1993), but also of the scales of structure in active galaxies (e.g. 
Saust 1994, Lewis et al. 1997). Statistical studies of numerical simulations have revealed that long term 
light curves of microlensed quasars will additionally probe the mass function of the microlensing objects 
(Lewis & Irwin 1996). 

In this paper, a complementary probe of these properties, via the identification of microlensing induced 
"centroid-shifting" in the images of macrolensed quasars, is presented. The paper begins with a description 
of gravitational microlensing due to stars in external galaxies (Section The microlensing model is 
described in detail in Section ||, with the results of this study being presented in Section ^. These are 
discussed in the context of space-borne interferometry in Section ^ and the conclusions of this study are 
presented in the final section. 



Since the discovery of the first gravitational lens system (Walsh et al. 1979), the action of the individual 
stellar mass objects on the light curve of the macrolensed object has be considered (Chang & Refsdal 1979, 
Young 1981). The problem introduces a natural scale-length for an isolated star of mass M, its Einstein 
radius 77, which is given by; 



where Djj are angular diameter distances between the observer (o), lens (1) and source (s) (Schneider et 
al. 1992). In this case the length is defined in the "source plane" at a distance Dqs from an observer. A 
point-like source passing within this projected radius of an isolated microlensing mass is amplified by a factor 
of at least 1.34, although the angular scale of multiple imaging induced by stellar mass lenses is ~ 10^^", 
far below present detection limits with optical telescopes. In practice, therefore, observations of microlensed 
systems are confined to the monitoring the image brightness fluctuations. The light curve induced by the 
passage of an isolated microlens possesses a very simple form, with a characteristic time-scale; 



2. Microlensing 



2.1. Background 




(1) 




(2) 
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where Vcs is the effective velocity of tfie source across tlie source plane (Kayser et al. 1986). In such simple 
cases, T can be used to characterize the mass of the microlenses (Wambsganss 1992). 

For a distant multiply-imaged quasar, whose light can shine through the inner regions of the 
macrolensing galaxy, the optical depth to microlcnsing can be substantial and the action of the individual 
lensing masses combine in a highly non-linear fashion, with the light curve of a background source exhibiting 
complex variability, including strongly asymmetric fluctuations (Paczyhski 1986, Wambsganss et al. 1990). 
In such a regime, the time-scale of individual "events" no longer reflects the time-scale given by Equation |^; 
such time-scales reflect the time taken for a caustic to sweep across the quasar continuum source and, as 
the details of the caustic network are due to the sum effect of the ensemble of microlenses, the properties of 
the lensing bodies must be determined from statistical considerations (eg Lewis & Irwin 1995,1996). 

Chang (1984) demonstrated that the ampliflcation of a source is dependent on its size relative to the 
Einstein radius of the characteristic lensing mass. At caustic crossings a point source is amplified by an 
infinite amount, but any physical extent leads to a finite amplification. Sources whose scale size is only a 
fraction of an Einstein radius can be amplified by large factors, while sources whose scale is much greater 
than an Einstein radius are amplified by a negligible amount. 



2.2. Q2237+0305 

Discovered during the CfA redshift survey (Huchra et al. 1985), the Q2237+0305 system consists of 
four images of a background quasar, z ~ 1.69, separated by < 2"0 (Yee 1988). The images are centered on 
the core of the lensing galaxy, a bright nearby, z ~ 0.04, barred spiral. The galaxy light in the inner regions 
of the Q2237-I-0305 system is dominated by stars in the central bulge of the galaxy (Yee 1988), and it is 
expected that optical depth to microlensing is substantial. This lensing galaxy is, however, unique among 
known gravitational lens systems in that it is located very close to the observer and thus the light -travel 
time differences to the four images are only ~ 1 day. Intrinsic variations in the luminosity of the quasar 
thus manifest themselves in all four images within ~ 1 day and differential variability between images 
with time-scales » Iday reflects the effects of microlensing. It is this effective decoupling of the observed 
time-scales for intrinsic variability and microlcnsing by stars in the galaxy that allowed photometric 
monitoring of Q2237-I-0305 to produce the flrst detection of a microlensing event (Irwin et al. 1989, 
Corrigan et al. 1991). Microlensing induced variability is a prominent feature of light curves of almost a 
decades monitoring of this system (0stensen et al. 1996). 
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However, for Q2237+0305, the Einstein radius (Equation [^) for a star of Solar mass, projected into the 

_ 1 

source plane, is 770 ^ 0.04ft,yg' parsec. (Through-out we adopt a standard Freidmann- Walker Universe, with 
51 = 1 of matter distributed smoothly. The cosmological constant. A, is assumed to be zero). Within the 
framework of the "standard model" for the central regions of quasars (Rces 1984) this microlensing scale is 
smaller than the extent of the broad emission line region, 0.1-1.0 parsec, but significantly in excess of the 
size of the ultraviolet-optical continuum-producing region, ^ 10^^ parsec. Thus, differential microlensing 
amplification of the quasar continuum and emission line producing regions by objects of stellar mass 
within the lensing galaxy of the Q2237+0305 system is expected. This is observable as a time-dependent 
variation in the relative strengths of the broad emission lines and the continuum which is coupled with 
microlensing-induced photometric variability (Lewis et al. 1997). 



2.3. "Improper Motions" 

Williams and Saha (1995) addressed the degree of centroid shifting during microlensing at substantial 
optical depth. For this study, a new model for mass distribution was constructed with non-parametric 
techniques, recovering both the relative image positions and magnitudes. Unlike previous models of this 
system, which predict a single image at each of the image positions (with each image consisting of a myriad 
of micro-images), the images in this model comprise of several "mini-images", separated by ^ 0"05. These 
mini-images suffer differing convergence and shear and, therefore, undergo differing amplifications. The 
close proximity of these mini-images implies that they are not resolvable with current optical telescopes 0, 
appearing as a single, point-like image whose central position is the centroid of the light distributions of 
the mini-images. 

As the light from each of these mini-images shines through the lensing galaxy it can be influenced by 
the action of microlensing bodies. This causes the brightness of each mini-image, and therefore the centroid 
of the light distribution, to vary over time. For the quadruple lensing system, Q2237-f0305, assuming the 
lensing galaxy possesses a transverse velocity of V = 600km s~^, observations on time-scales of one year 
have a 50% chance of revealing shifts of <L 0"003, while on ten year time-scales there is similar chance 

^Such image splittings may be accessible through radio observations with the VLBI, which are capable of 
probing milli-arcsec scales, although Q2237-I-0305 is not a powerful radio source with observed flux densities 
of 40 - 80 ^ Jy at 3.6 cm (Falco et al. 1996). 
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of ^ 0"01, and a 15% chance of ^ 0"03. Centroiding of images to a fraction of a resolution element, 
typically ^ 0"01, is now routine and Williams and Saha (1995) propose that observed centroid shifts of 
these magnitudes, observed in Q2237+0305, could be due to their proposed mechanism. 

It should be noted, however, that the splitting of a macro-image into several mini-images which 
experience differing microlensing characteristics is particular to the modeling procedure of Williams and 
Saha (1995); other techniques predict a single macroimagc comprised of micro-images which have a typical 
separation of micro-arcsec. In this paper we examine this regime by centroid shifts introduced by the 
variation in brightness of these micro-images as stellar field passes in front of the macrolensed source. 

3. Microlensing Model 
3.1. Method 

A number of techniques have been developed to numerically study the action of microlensing objects 
on distant sources (Young 1981, Paczyhski 1986, Witt 1993, Lewis et al. 1993). The "work-horse" of 
numerical simulation has been the backwards ray-tracing method (Kayser et al. 1986, Wambsganss et al. 
1990) in which a regular grid of rays is fired from an observer into a field of microlensing masses lying in the 
"image plane" . The gravitational lens action of these masses act to deflect the ray. The ray is then "traced" 
until it impacts the source plane. Although the input rays are regularly distributed, the non-linear action 
of the microlenses ensures the distribution of rays over the source plane is not (Figure |l|) and the density of 
rays at a point in the source plane, relative to the density of rays over the image plane, is a measure of the 
amplification a source at that position would experience. 

To achieve sufficient resolution over the source plane, however, a large number of rays needs to be 
traced through the field of microlensing masses. Compounding this, the direct calculation of the deflection 
angle suffered by a ray due to a large number of these masses is computationally expensive and simulations 
involving direct calculations would take a prohibitive time to complete. To alleviate this, the gravitational 
microlensing "force" is approximated with the use of tree codes (Barnes and Hut 1986), allowing realistic 
simulations to be performed on desk-top computers. 

Such a code is employed in this study. The source region was chosen to be 2.5 Einstein Radii for a Solar 
mass, and pixelated to 256^. With this, the area of source plane under consideration is ^ 3.2 x lO^^h^g'cm, 
and each source plane pixels is 1.2 x lO^^hy^'cm along a side. A ray fired through the lens plane at 
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(xi,yi) will experience a deflection due to the microlcnsing masses and will intercept the source plane at 
a position (xs,ys)- Considering a single source, at position (xR,yR), it can be seen that once a grid of 
Nr rays has been traced through the lensing plane, the centroid of the macroimage, due to the individual 
microimages, is at; 



with a similar equation for the y co-ordinate. Here, the function p describes the surface brightness 
distribution of the source. For the study presented here, a source is assumed to be a disk of radius R, with 
a uniform surface brightness distribution. More general source profiles can be employed, although for this 
study the scale of any centroid shift is of interest and this will not depend heavily on the detail of the source 
profile. 



As with the earlier work of Williams and Saha (1995), we consider image centroiding in the multiply 
imaged quasar, Q2237+0305. Since its discovery, a number of models, both parametric and non-parametric, 
for the lensing mass distribution in this system have been developed. Using ground-based (Kent and Falco 
1988, Schneider et al 1988), and later, Hubble space telescope data (Crane et al 1991, Rix et al 1992, 
Wambsganss and Paczyiiski 1994, Schmidt et al 1998), these models can, to within the observational error, 
reproduce the image configuration. Many models, however, fail to accurately reproduce the observed optical 
flux ratios, although there is broader agreement with the ratios from the larger radio emitting region (Falco 
et al. 1997). It has been suggested that this discrepancy is due to the action of differential microlcnsing, 
with the smaller UV-optical emitting continuum source being more susceptible to the induced fluctuations 
(Schmidt et al 1998), although the effects of differential reddening between the images has also been 
considered (Nadeau et al. 1991). 

Although these models vary greatly in the details, those possessing parameterized, elliptical forms for 
the mass distribution, without the need for an additional external source of shear, tend to predict similar 
optical depth and shear parameters at the positions of the images ^ For simplicity we employ the "best-flt" 
model of Rix et al. (1992), which employs an elliptical r^i mass distribution with an unresolved core. The 

^It should be noted that this is also the case for the brightest of the three "mini-images" that comprise 
image B of the non-parametric model of Williams and Saha (1995). 
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3.2. 



Model Parameters 
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microlensing parameters at the positions of the images are summarized in Table 1 



Image 


(7* 


7 




A 


0.41 


0.47 


7.86 


B 


0.38 


0.43 


5.01 


C 


0.65 


0.68 


2.94 


D 


0.59 


0.56 


6.87 



Table 1. The macrolensing parameters, the surface density of matter, (T,, and large scale shear, 7, at the 
positions of the quasar images in Q2237+0305. The values are the "best" model parameters of Rix et al. 
(1992), and Hth is the "macrolensing amplification". The image designations are those of Yee (1988). 

For each image, two simulations were performed. Firstly, the microlensing stars were distributed with 
a uniform mass function, whereby all objects have a mass of IMq, and secondly according to a Salpeter 
mass function of the form 

p{m) dm oc m'~^'^^ dm (4) 

between the limits O.OSMq < M < IOMq, with (M) = 0.163Mq. 

Three source radii were employed, R = 5 x 10^*, 5 x 10^^ & 5 x lO-'^^hyg'cm which are of order the 
scales of structure expected for continuum emission from an accretion disk at a quasars core (Ress 1984). 
For these sources, a total of 1275 rays per pixel were traced through the image plane, corresponding to 
a total number of rays, collected in the source plane, of 84 million per simulation. To ensure that the 
vast majority of micro-images are considered, the total number of rays traced per simulation was several 
times greater than this value. Considering the smallest source, and assuming a maximum deamplification 
of 2 magnitudes, this number of rays ensure a maximum error of less than 10% in the derived image 
amplification. 

■^There is an error in the equation for the shear parameter given in Table 3 of Rix et al. (1992). This 
should read 7 = y^(l — t)^ — 1/Ampl, where t is the microlensing optical depth (= cr*) and Ampl is the 
image amplification. Taking the Rix et al. (1992) values for the optical depth and amplifications, the shear 
values at the positions of the images should be 7^ — 0.43, 75 = 0.43, 7c — 0.69 and jd = 0.62. The 
differences between these and the employed microlensing parameters do not effect the overall conclusions of 
this paper. 
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4. Results 

Figures || and ^ present the results of this study for the constant Solar mass distribution and Salpeter 
distribution respectively. Each figure consists of four main panels, one for each image in Q2237+0305. Each 
panel presents, on the right and left side, two trajectories across the amplification map, each perpendicular 
to one another. With this, one trajectory is aligned with the dominant component of the shear. The 
statistical properties of induced variability a source experiences as it travels behind a microlensing region 
depends on its orientation with respect to this shear component, with more rapid variability occurring 
as the source moves along the dominant shear component, although the small scale of the source plane 
under consideration does not make this direction immediately apparent. Each Figure presents data for 
the three source sizes, the black line representing the smallest source, the light grey line, the largest. The 
intermediate source is represented by a grey line, although it follows the variability of the smallest source 
to a high degree, with significant differences being apparent only in the vicinity of caustics. 

The top figures in each panel present the image light curves, showing, in magnitudes, the fluctuations 
about the "macrolensing amplification" (Table |^). As expected, the smallest source undergoes the most 
rapid and complex variability, while the largest source exhibits long term trends. Below each light curve, 
two frames are presented. The first presents the degree of centroid shift in the direction perpendicular to 
the shear. It should be noted that this shift is expressed in units of 15 angular Einstein radii. In all the 
frames this component is negligible on this scale, with several showing low, long term trends. The bottom 
frames show the centroid shift in the direction of the major shear component, and these display considerable 
fluctuations on both short and long time scales. As with the image light curves, these fluctuations arc more 
severe for the smallest sources, with the larger source following the large scale trends in the variability. 
Any large change fluctuation in the centroid position is correlated with a fluctuation in the total image 
brightness, although the converse is not true and all fluctuations in the image light curves do not result 
in large centroid shifts. This indicates that large changes in the macro-image centroid are due to the 
production or disappearance of bright images during the crossing of a caustic, although some crossings 
produce image pairs close enough to a current image centroid to produce no shift. 

To illustrate this effect in more detail we consider a small displacement of the quasar source, and 
examine the resultant change in the image configuration. Utilizing the simulation for Image A with all stars 
possessing a mass of 1 Mq the source was placed at two points straddling a caustic. These source positions 
are illustrated, superimposed upon the caustic network, in Figure |^. The position of the source, relative 
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to this caustic network, changes by only by a smaU amount (AEg = 0.2ER), but as seen in the upper-left 
panel of Figure ||, the resultant change in the image position is substantial, being AEj ^ 3ER for the 
largest source, and greater for the smaller sources, with AEj ^ lOER. The distribution of images at these 
two source positions are presented in Figure ^ The 'low' resolution of this Figure allows us to consider in 
detail the effects on the largest source only. At each of the time steps, the distribution of micro-images is 
dominated by three large (bright) images, which lie along the direction of the dominant shear component. 
If all the matter in the galaxy were smoothly distributed the resultant image centroid would be on the left 
of the frame (denoted by the symbol). The smaller to larger circles also presented in the plot denote 
the centroids of the smallest to largest source size under consideration. Between the frames, the left-hand 
images fade, as the source moves from the region of high amplification, through the caustic, and into the 
deamplification region. One of these images (far-left) appears to be the pair of images merging due to the 
caustic crossing. Due to its extent, the wing of the largest source still overlaps the caustic at the latter 
time step and possess an appreciable image on the far left of the frame. The smaller sources, which would 
completely crossed the caustic by this point, have lost these bright images. 

The large-scale distribution of micro-images is also apparent in Figure 3 of Wambsganss (1992). This 
presents the area of the image plane from which rays are traced into a particular square region of the source 
plane, as compared to a similar region which would map directly onto the source plane area if the material 
in the microlensing stars was smoothly distributed. The area of the image plane is effectively the image of 
the square source plane region and is composed of five major components, although, due to the small region 
of this image plane under consideration, one of images extends beyond the figure. It is apparent that the 
centroid of these images will not coincide with that of the image created by a smooth matter distribution 
by a few Einstein radii. This frame represents a single snap-shot of possible micro-image configurations 
and more extreme fluctuations are likely. 



Although not indicative of the time scales of individual events, the crossing time of an Einstein radius 
(Equation ||) still provides a temporal measure of microlensing activity. Following Kayser et al. (1986), this 
crossing time for a Solar mass star in Q2237+0305 is 



4.1. Time Scales and Physical Shifts 



T 




V600 



8.7 



hrs'yrs , 



(5) 
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where zi is the redshift of the lensing galaxy, and vi — GOOvgookni s^^ is its departure velocity from the 
Hubble flow, perpendicular to the line-of-sight. Figures || and |3| reveal that the smaller sources experience 
centroid shifts on much shorter time scales than this, with substantial shifts occurring on <yr periods, 
superimposed on longer term variability. 

At a redshift of z = 0.04, the lensing galaxy in Q2237+0305 lies at a distance of 147hy5Mpc, and the 

1 c - 

angle subtended by the Einstein radius of a Solar mass star (0.0052h7g^pc) is 7.3 x 10^ h^g". This implies 
that the unit of centroid shift in Figures ^ and || is 110 x lO^^hyg". 

Considering current limits of optical astrometry, such centroid shifts appear to be well below detection. 
In the next Section, however, we consider the astrometric sensitivity of a forthcoming space craft, the Space 
Interferometry mission. 



5. Space-borne Interferometry 

Over its history, much of astronomy has focused upon the measurement of distances to objects in 
the universe. Other than a small number of "direct" methods, such as gravitational lens time-delay 
measurements (eg Kundic et al. 1997, Schechter et al. 1997) or the detection of a cluster's Sunyaev- 
Zeldovich decrement (eg Myers et al. 1997), most rely on a "ladder" of techniques to determine the 
distance to extra-Solar objects. The lowest rung of this ladder requires the determination of distance to 
nearby objects, allowing a calibration of their physical properties, and this can be achieved with a purely 
geometrical measure, the parallax. This method does, however, require accurate determination of the 
"shift" of objects with respect to the extra-galactic background as the Earth orbits the sun. The effects of 
atmospheric seeing limits such measures, allowing a parallax determination of fractions of an arcsecond. 

Astrometric measurements will be revolutionized with the development of space-borne optical 
interferometric techniques (Loiseau and Malbet 1995). Due to fly in 2005, the Space Interferometry Mission 
(SIM) promises to surpass all previous space-borne astrometric missions [e.g. Hipparcos (Ferryman et al. 
1992)] by several orders of magnitude. Its current operational parameters will provide centroiding of a 20**^ 
magnitude stellar image to 4 x 10^^" with 14 hours of integration Such accuracy will provide parallax 
determinations of Cepheids out to ~ 5kpc, and will allow the detection of 'wobbles' in stellar proper 

"^For a more comprehensive review of the Space Interferometry Mission see 



http : / /huey ■ jpl ■ nasa . gov/sim/ 
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motions due the presence of Earth mass planets about nearby stars. Recently, Miyamoto and Yoshii (1995) 
proposed to use SIM to observe the 20/i" shifts expected during a microlensing event in the galactic halo, 
providing more constraints on the lensing geometry. As halo microlensing events are typically due to an 
isloated MACHO, astrometric observations with SIM will determine the mass, distance and proper motion 
of the microlensing body, even if the MACHO is too faint to be detected directly (Paczynski 1998). 

It is immediately apparent that the scale of microlensed induced image shifts in the centroids of 
macrolensed quasars {niR ~ 17 — 18.5) are of sufficient magnitude and occur on short enough time scales to 
place then are well within the operational capabilities of SIM. As with photometric monitoring, however, the 
complex nature of the caustic distribution at high optical depths will prevent the association of a particular 
image shifting "event" with an individual mass (eg. Wambsganss 1992); rather the statistical properties 
of a number of image shifts, such as their duration and degree, will provide measure of a characteristic 
lensing mass and scale-size for the continuum emitting region in the core of the quasar. Importantly, the 
dominant direction of image shifting is coincident with the direction of the dominant shear component, and 
the determination of any shifts perpendicular to this direction will provide a measure of the relative values 
and orientation of the two amplification eigenvalues. A non-detection of any perpendicular shifts still 
reveals the orientation of the amplification eigenvalues, but provides a limit on their relative values. The 
identification of these components will provide powerful additional constraints on any macrolensing mass 
model for this system, leading to a more robust measure of the microlensing optical depth and shear at the 
position of the images. 



6. Conclusions 

This paper has investigated the role of microlensing in inducing image shifts in the centroids of 
the star-like images of multiply imaged quasars. It was found that for the quadruple quasar system, 
Q2237+0305, that substantial shift of ~ lOO/x" are possible on time scales of months. Although not 
observable with current telescopes, such shifts will be detectable in the next decade with the proposed 
Space Interferometry Mission. This initial study has revealed that the characteristics of centroid shifting 
are sensitive to the size of the source undergoing microlensing, providing an independent probe of the scales 
of structure in the continuum emitting region of quasars. Similarly, these characteristics should depend on 
the mass function of the microlensing objects, although a proper determination of any influence will require 
statistical study of an ensemble of simulations. 
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We wish to emphasize, however, that the Ught curves and image shifts presented here are drawn from 
only single realizations of the microlensing in the images of Q2237+0305. The results must, therefore, 
be taken as purely illustrative of the effect. Currently, an ensemble of simulations are underway; once 
complete, these will be employed in a study of the statistical properties of image centroid shifting and their 
dependence upon the scale of the source and mass function of the microlensing bodies. 
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Fig. 1. — An amplification map for Image A in Q2237+0305 (cr* =0.41,7 = 0.47). The grey-scale is 
a measure of (de) amplification and is in units of magnitudes. Note that a AM is with-respect-to the 
brightness the macro-image would appear if all the matter was smoothly distributed. The map has an 
extent of 2.5 Einstein radii for a Solar mass star along each axis, and the microlensing stars were distributed 
with a uniform mass of IM©. The circles represent two positions of the source which are discussed in more 
detail in Section ^. 

Fig. 2. — The results for the individual images of Q2237+0305. All the microlensing stars possessed a 
mass of 1M0. For each of the images, (A,B,C,D), two sets of panels are presented, each representing source 
trajectories that are either aligned or perpendicular to the direction of the dominant shear component. Each 
plot is a function of the source position along this trajectory in units of Einstein radii for a Solar mass star 
in the source plane. For each orientation, three plots are given. The top panels for each image presented the 
light curve (in magnitudes) for the smallest source (darkest line) through to the largest (lightest line). The 
middle panel presents the centroid shift in the direction perpendicular to the major shear axes, while the 
lower panels are the centroid shift in the direction of the major shear component. Note that in these lower 
plots that shift is expressed in units of 15 angular Einstein radii for a Solar mass star. 

Fig. 3. — As for Figure H, but for the stars distributed with a Salpeter mass function. All the axes are in 
the same units as Figure ^ including the 'x-axis' which represents angular Einstein radii of a Solar mass 
star. 

Fig. 4. — The distribution of micro-images for the two source positions presented in Figure |l|. The units 
along the axes are angular Einstein radii of a Solar mass star. The dark grey areas are the images for 
the largest source (R = 5 x lO^^h^g'cm), while the lighter grey area represents the intermediate source 
(R = 5 X lO^^hyg^cm). Note that the resolution of this image does not reveal the full image structure for this 
smaller source, and for the R = 5 x lO^^hyg'cm source. The symbol indicates the position of the center of 
the image if all the lensing matter was smoothly distributed. Each circle corresponds to the centroid of the 
micro-image brightness distribution for each source, the smallest circle being for the smallest source etc. 
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